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Based on an exact expression for the self-energy of the Jahn- Teller polaron, we find that symme- 
try of pseudospin rotation makes the vertex correction much less effective than that for the Holstein 
polaron. This ineffectiveness brings about a smaller effective mass m* and a quantitatively differ- 
enent large-to-small polaron crossover, as examined by exact diagonalization in a two-site system. 
In the strong-coupling and antiadiabatic region, a rigorous analytic expression is found for m*. 
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It is well recognized that both the double-exchange 
mechanism and the strong electron-phonon interac- 
tion, specifically the Jahn- Teller (JT) effect on doubly- 
degenerate Cg orbitals coupled with two degenerate vi- 
brations (the E ® e case) at each Mn'^+ site, are essen- 
tial ingredients to bring about the colossal magnetoresis- 
tance (CMR) in manganese-oxide perovskitesJja Thus 
the theories oruGMR need to include these ingredients 
simultaneouslyJ3~l3 This complicated situation compels 
some theories to neglect kinetic energies of ions and oth- 
ers to treat the JjT polaron in a way similar to the conven- 
tional polaron.EfLI In either way, characteristic features of 
the JT polaron do not emerge from those theories. 

In fact, in spite of a broad interest in its role in su- 
perconductivity, studies on the JT effect in an itinerant 
electron system are limited, probably because the word 
"the JT effect" often implies strong lattice deformations 
and- a localized electron associated with them. Hock et 
aZ.El considered the simplest case, namely, the E® (3 JT 
polaron which, unfortunately, possesses a too simple in- 
ternal structure to provide qualitatively different features 
from those of the Holstein polaron.El The second simplest 
E®e case was treated by, Fabrizio and TosattiO as well as 
Benedetti and Zeyher,E3 but both works addressed only 
localization in the strong-coupling region. 

In this paper, we provide a new aspect to this problem 
by making a comparative study of the E ® e JT polaron 
with the Holstein polaron based on the knowledgej-at- 
tained after forty-year's investigation into the latter .Ej 

We have obtained the following results for the system 
specified by the two parameters, t = t/ujQ and a = i?jT/t^Oj 
where t, wq, and -Ejt are the energies corresponding 
to bare electron transfer, bare phonon, and Jahn- Teller 
stabilization, respectively. (1) Based on an expression 
for the self-enerffii_derived by a similar method for the 
Frohlich polaron,Efl we find that the vertex correction is 
much less effective in the JT polaron than that in the 
conventional polarons due to aJocal conservation law im- 
posed on the JT HamiltonianJlj (2) This ineffectiveness 
leads us to a smaller effective mass, as shown by an ex- 
plicit expression for the JT polaron-mass enhancement 
factor as yj2/'Ka e" in the antiadiabatic {i <C 1) and 
strong-coupling (a ^ 1) region. (3) The large-to-small 
polaron crossover is examined by exact diagonalization 



(ED) in a two-site system on the ground that the ED 
calculation on small clusters is very effective for a':^!. 
We find that the crossover occurs at i « a~\foL, indicating 
that JT polarons are more mobile than Holstein ones. 

Let us start with a single_p0e center at site j described 
by the Hamiltonian iJj adlj 

-A[(7ja(4„dj6 + 4dja)+gj6(4,rfja-4rfj6)] + ffjv, (1) 

with A — u}Q\fl^EyY, where dja and djf, represent electron 
annihilation operators for the two degenerate orbitals, 
(7ja and gjf, are the two local JT distortions, and i/jv 
is the harmonic Hamiltonian for the vibrational modes. 
In polar coordinates as gja = (7j sin6'j and gjf, = gj cos^j, 
the energy eigenfunction for Hy^(^q-fi-^nV) ^ satisfying 
}l-^^\nl)—ujQ(n + \)\nl) is given byO 

(gj^jln/) = iV,,pi^(-p, |;| + l,zj)zf/2e-^j/2e''«^ (2) 

with iV,,p = (_i)(l'H)/2^^^(|;|+p)!/^p!/|^|!^ = ^^^2^ 
and n—\l\-\-2'p, where F{—p, Zj) is the confluent hy- 
pergeometric function, I is an integer, and p — 0, 1, 2, • • •. 

In terms of boson operators, aj and bj, to represent q^a 
and Qjb in second quantization, Hj is rewritten as 

Hj = wo\/2^[(aj+-6j)Cj\cjx + (aj-fej+)c+Cjt] 

+ujo{afai + b+bi + l), (3) 

where pseudospin index a{= ±1 =t or J,) for elec- 
tron operators is introduced through the relation cjo- = 
{dja+i<jdjb)/ V^- Note that second-quantized representa- 
tion for phonons is not unique due to SU(2) symmetry 
in ifjv We have chosen it in such a way as to diagonalize 
both _Hjv and lj = ~id/d6y Then we obain \nl) as 

^ n + l n-i 

\nl)~ — ^ hf ^ I vacuum). (4) 

V[(n+0/2]![(n-0/2]! ^ ^ ' 

This phonon representation is a key step to obtain ana- 
lytic expressions in Eqs. (||) and (|l^) as well as a clear 
view of less effectiveness of the vertex correction. 

Because of the symmetry of pseudospin rotation, the 
operator Lj, defined by Lj = Oj^aj— 6j^5j— (Cj^Cj|— Cj^Cj|)/2, 
is conserved as easily checked by [iJj, Lj] =0. 
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For a single electron at site j, eigenvalues for Lj are 
half-integers and each energy level is doubly degener- 
ate corresponding to ±|Lj|. In general we can give the 

ground-state wavefunction "^^^^ only numerically, but for 
large a we find an analytic expression as 



(0) 



(5) 



for Lj ~ l/2,ll3 where Jo{x) is the Bessel function and 
Ii{x) its modified form. The corresponding energy Eq is 
given as « (—a -|- 1/2 -I- l/16a)ajo ~ — i?jT- 

Now we consider a lattice composed of N JT centers 
for which the Hamiltonian is given as i?jT = Ht + J^j j 
where Ht describes the transfer energies between nearest- 
neighbor JT centers as 



Ht 



^ E E tff'idfdyy + 4^,dj^). (6) 



<jj'> 7=a 7'=a 



For simplicity, we take t^-y' =5^~fit in the following. Then 
Eq. (o) can be rewritten as 



<jj'> 



'jt) = E'^kCi^^Ck^, (7) 



kcr 



where c\^ct{—N~^^'^ e~*-' *'cjo-) is the Fourier transform 
of Cj(j and £k represents its bare dispersion relation. Note 
that the operator L defined by i = ij is conserved, 
namely, [-ffjT, L]—Q in this choice of t^^i. 

The thermal one-electron Green's function G\^a{iuird 
with u)n a fermion Matsubara frequency is defined byllj 



Gkcr(jw„) 



/ dTe^'^"^Gk.(T), 
Jo 



(8) 



with (3 = T-i and Gkirlr) = -(TVck^(T)ci^^). We first 
consider 9Gk<j (t) / dr to derive an equation of motion 
which relates Gko-lT) with the electron-phonon corre- 
lation function (?V Z]q{[a^('^)-^-q('^)]ck+q-<T(T)c+^}). 
Next we derive a similar equation of motion for this 
correlation function in order to eliminate phonon op- 
erators in the expressions other than the bare phonon 
Green's function which is the same for both phonons as 
i?o(*i^m) = 2u;o/[(ia;m)^ — wo^] with tOm a boson Matsub- 
ara frequency. Then we arrive at an exact expression for 
the self-energy I]kcr(*'^n) as 

cj„/ k' 

xGk'-cr(it^n')A-o-o-(k', iw„'; k, zw„), (9) 

where the vertex function ACT'(T(k', ia;„'; k, icjn), a key 
quantity in this expression, is found to be 

Gko-(«^ri)Gk'o-'(«a;„')Acr'o-(k',iw„';k,2Cj„)= / dre'""'^ 

Jo 

X /''dr'e'^— "')-'(T.Ckv'(T)5£:^k(T')4.>, (10) 
Jo 



with <S'£/_5c = Z^k" '^k"-*'-k(T''^k"CT, reflecting the spinor na- 
ture of the problem. Equation (|9|) serves as a firm basis to 
study the JT polaron in the Green's function approach. 

Quite an analogous-j-esult has been obtained for the 
conventional polaronEJ For the Holstein model specified 
by the Hamiltonian iJjj as 



i^H = E ^kCi^^Cko-+Wo\/2^EK^+°j)^jV^j° 



k(T jo- 

+a>oE(4«J + l/2)' 
j 



(11) 



where a in this case refers to "real" spin index, Ek(j(jw„) 
is given in the form of Eq. in which Gk'-cr(*^n') 
and A_o-o-(k', ia;„'; k, zwn) are, respectively, changed into 
Gkvli^n') and Ac(k', io-in'; k, ia;„) the charge vertex 
function, defined through Eq. (|l^) with 5'^/,^^. replaced 

by the charge operator pk'-k = Sk'V" ct'-*'-ka"'^k"o-" 
due to the scalar nature of the Holstein system. 

The diagram to represent Eq. (|^) is given in Fig. 1(a), 
in which we introduce the vertex F by eliminating im- 
proper diagrams from the vertex A. The expansion se- 
ries for F in terms of a is shown in Fig. 1(b). If we as- 
sume that G-k_ij{iuJn) is independent of a and emply the 
Migdal's approximationE2l in which only Fq is retained 
for F, namely, F = 1, there exists no difference in the 
self-energy between JT and Holstein systems. 




^2d J^2e " ''Z 

FIG. 1. (a) Self-energy in diagrammatic representation. 
Thick solid and thin dashed lines indicate, respectively, the 
electron Green's function and the bare phonon propagator, 
(b) Expansion series for the vertex F up to second order in a. 

There is, however, an important difference in the ver- 
tex correction. In contrast to the Holstein system, the 
corrections represented by the diagrams Fi, F2a, ■ • ■ , F2e 
is seen to vanish in the JT system by merely considering 
the pseudospin assignment together with the direction of 
phonon propagators, because Eq. (||) dictates that the 
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JT-phonon exchange interaction works only in the pseu- 
dospin exchange process between electrons with opposite 
pseudospins. Physically both electrons and phonons in 
the JT system are associated with a notion of clockwise- 
or counterclockwise- "rotation" around each JT center 
and electrons interact with phonons only when the to- 
tal rotation is conserved. In this sense, the vanishment 
of these vertex corrections is due to the local-rotation 
conservation law. This law allows only processes such 
as the one represented by T2f for T. Similarly, all the 
third-order vertex corrections vanish. Ineffectiveness of 
the vertex correction widens the applicable range in a 
of the Migdal's approximation in the JT system and it 
leads us to the smaller polaron mass enhancement fac- 
tor m*/m than that in the Holstcin model in which the 
correction Fi is known to enhances m* /m as a increases. 

The above perturbative approach is not useful in dis- 
cussing a small polaron or polaron localization in a^ite. 
According to the studies on the Holstein model,ll3 an 
ED calculation in a two-site system provides qualitatively 
correct and quantitatively fair results for the small po- 
laron in the strong-coupling region (a > 1), irrespective of 
the value of t. Thus we shall make a similar analysis of a 
single electron in the JT model with iV = 2 in which the 
eigenvalues of the conserved quantity L are half-integers 
and each energy level is doubly degenerate. 

Let us consider the antiadiabatic region (i <C 1) first. 
At a^l, both the ground and first-excited states belong 
to the sector of \L\ = 1/2. Using in Eq. (||), their 

wavefunctions are written as ^'j. « (^'^"•' ± ^!^'^)/\^ 
for L = 1/2 with the corresponding energies E± = 
Eq ± t/[lQ{a) + Ii{a)]. The energy difference, E+ — E^, 
can be used to estimate the polaron bandwidth in a crys- 
tal and thus its ratio with the bare value, 2t, determines 
the polaron effective mass through the relation 
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E. 



E^ 



2t 



Io{a)+Ii{a) 



(12) 



This result should be compared with e ^" the Holstein's 
famous resultti for the system defined in Eq. ( pi] ) . 

We resort to ED calculations to obtain m/m* through 
the numerical evaluation of E± for arbitrary a. The con- 
servation of L helps reduce the number of expansion bases 
for phonons considerably. We plot the calculated m/m* 
for both JT and Holstein models in Fig. 2 in which t is 
taken as 0.2, although the result itself does not depend on 
t provided that it is much smaller than unity. (The result 
for the Holstein model hardly changes from the analytic 
result e~^" in the whole region of a.) For small a, both 
models give essentially the same m/m* as implied by the 
previous weak-coupling analysis. For large a, however, 
there is a difference in m* /m which is more than orders 
of magnitude for a> 1, indicating that the JT polaron is 
quite mobile compared to the Holstein polaron. 
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FIG. 2. Polaron mass reduction factor m/m* for the JT 
(the solid curve) and the Holstein (the dashed curve) models 
with each analytic expression in the strong-coupling region. 

Next we make a semiclassical argument on the adia- 
batic region (t^l) by considering the adiabatic poten- 
tial [/ad for given phonon variables, {(/i^i; (72^2}- Since it 
was calculated previously in connection with the Berry 
phasejH^ we just give the result here as 



'2)]+c?^o%i'^q2'fY"}"\ (13) 



X C0S(fi — f 



with q^ = qi^ + q2^. This potential has rather simple fea- 
tures; if the adiabaticity parameter \=a/t — E}^ /t is less 
than unity, C/ad has only one minimum in {91^1; 92^2}- 
coordinate space, implying no symptom for a small po- 
laron. On the other hand, it is a double- well potential for 
A>1 with the energy barrier A = (q:wo/2)(1 — A^^)'^. If 
the largest zero-point energy of phonons Azdo (which is 
loq/2 in this case) is smaller than A, localization leading 
to a small polaron occurs. Thus the condition A ^ Azcro 
provides the criterion to obtain a small polaron as 



t ^ a ~ ^/a, with a > 1. 



(14) 



A similar argument has been done for the Holstcin model 
described in Eq. (|l^) for which a double-well potential 
appears only when A>l/2 with A = acjo(l^l/2A)^ and 
Azcro = (^0 /2) \/l — 1/ 4A^1, leading to the criterion 
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(-2^) 



2Ay ^ 2 



1 - 



4A2 



(15) 



This condition cannot be reduced to such a simple form 
as that in (^^, but clearly it is much less restrictive than 
( p^ ) for the small-polaron formation. 

Finally we make a more quantitative argument on 
the large-to-small polaron crossover based on the exact 
ground-state wavefunction obtained by the ED calcu- 
lation. We evaluate two quantitites, "the transfer ampli- 
tude per bond" T= (I-qI Ea(C]*aC2^ + 4aCi<T)|*o>/A^bond 
with the number of the bond A'^bond = 1 and "the inter- 
action amplitude per site" /= J2j [{'^^ ~bj)c'^^Cji + 
(flj -b'^)c'^iCj^] \^q)/N with A = 2. Then we measure 
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"itineracy" by the ratio |T//|, because the ratio must be 
large for an itinerant polaron. 

Contour pfots for \T/I\ in (f, Q;)-plane are given in 
Fig. 3. The result for the Holstein polaron indicates 
that the semiclassical criterion for the small-polaron for- 
mation corresponds to the condition \T/I\ « 0.5. More 
or less the same result is obtained for the JT polaron 
for which Eq. (^) is well represented by the condi- 
tion \T/I\ K, 0.6. In either way, we can conclude that 
the large-to-small polaron crossover occurs at around 
|r//|«0.5 — 0.6 and that a small polaron is much harder 
to realize in the JT system than the Holstein one. 
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FIG. 3. Contour plot for \T/I\ for (a) JT and (b) Holstein 
polarons. (Only the curves in the range 0.1 — 1.0 are shown to 
avoid too many curves.) The thick dotted curves correspond 
to the semiclassical criteria to divide large and small polarons, 
Eqs. Q and (|l|). 

Three comments are in order: (1) In the manganese 
oxides, the parameters are estimated as i « 0.2eV, wq ~ 
O.OSeV, and Ej^ « 0.2 - 0.4eV, leading to 2.5 and 
a ~ 2.5 — 5, which covers the crossover region according 
to Fig. 3(a). This is convenient to explain the observed 
CMR behavior. (2) The very large to* in the Holstein 
model is unfavorable pfor the bipolaron scenario for high- 
Tc superconductivityE3 In this respect, a smaller m* was 
suggested for the Frohlich polaron.O The same may be 
claimed for the JT polaron. (3) The electron-phonon 
coupling constant in H^i [Eq. (O)] is so determined as to 
give the same polaron effect as the JT case in the weak- 
coupling region for the proper comparison of vertex cor- 
rections. In this choice, the ground-state energy for i?H 



at i = is given as {—2a+\/2)u}Q which is about — 2i?jT 
at 1. Thus, if we make an alternative choice of the 
coupling constant as to give the same polaron stabiliza- 
tion energy in the strong-coupling limit, the difference 
in m* between the JT and Holstein models looks to be 
much reduced, but even in this choice, the JT polaron 
has smaller to* at least by the factor of Xj^fa. 

In conclusion, we have compared the E e JT po- 
laron with the Holstein one by using various theoretical 
techniques. Features of these polarons are exactly the 
same in the weak-coupling region, but they are different 
quantitatively in other regions due to the symmetry of 
pseudospin rotation; the JT polaron is more mobile than 
the Holstein one. 
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